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14. Abstract/Notes 


Many applications using LANDSAT images face a dilemma: the 
user needs a certain scene (for example, a flooded region), but that 
particular image may present interference or noise in form of horizontal 
stripes. During automatic analysis, this interference or noise may cause 
false readings of the region of interest. In order to minimize this 
interference or noise, many solutions are used, for instance, that of 
using the average (sunple or weighted) values of the neighbouring vertical 
points. In case of high interference (more than one adjacent line lost) the 
method of averages may not suit the desired purpose. The solution proposed 
here is to use a spline-like algorithm (weighted splines). This type of 
interpolation is simple to be computer implemented, fast, uses only four 
points in each interval, and eliminates the necessity of solving a linear 
equation system. In the normal mode of operation, the first and second 
derivatives of the solution function arc contini 'S and determined by data 
points, as in cubic splines. It is possible, however, to impose the values 
of the first derivatives, in order to account for sharp boundaries, without 
increasing the computational effort. Some examples using the proposed 
method are also shown. 
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NOISE CORRECTION ON LANDSAT IMAGES USING A SPLINE-LIKE ALGORITHM 


Luiz Alberto Vieira Dias (1) 
Nandamudi Lankalapalli Vijaykumar (2) 


ABSTRACT 

Many applications using LANDSAT images face a dilema: the user needs a certain 
scene (for example, a flooded region), but- that particular image may present 
interference or noise in form of horizontal stripes. During automatic 
analysis, this interference or noise may cause false readings of the region of 
interest. In order to minimize this interference or noise, many solutions are 
used, for instance, that of using the average (simple or weighted) values of 
the neighbouring vertical points. In case of high interference (more than one 
adjacent line lost) the method of averages may not suit the desired purpose. 
Immediate solutions that nay be thouoht of are polynomial or spline 
i' terpolations. The polynomial interpolation has the disadvantage of 
presenting oscillations (depending on the polynomial degree), while, in the 
second method, one has to solve a linear system that may be large depending 
on the number of points and consequently takes more time. The solution 
proposed here is to use a spline-like algorithm (weighted splines). This type 
of interpolation is simple to be computer implemented, fast, uses only four 
points in each interval, and eliminates the necessity of solving a linear 
equation system. In the normal mode of operation, the first and second 
derivatives of the solution function are continuous and determined by the data 
points, as in cubic splines. It is possible, however, tc impose the values of 
the first derivatives, in order to account for sharp boundaries, without 
increasing the computational effort. Some examples using the proposed rethod 
are also shown. 


1. INTRODUCTION 

In order to minimize the noise or interference ’n form of horizontal stripes in LANDSAT 
imagery, many solutions are used, like substitution. of the stripe pixel values for the average 
(simple or weighted) values of the neighbouring vertical points. This solution may not suit if more 
than one adjacent line is lost or if the pixel value variation is not linear (as observed, the 
variations are seldom linear). 

Sometimes the user needs a certain scene, and the particular images are noisy or with 
interference. It is interesting to have a nonlinear interpolator, that has not the disadvantages of 
the high degree polynomial interpolations (oscillations) or splines (need of solving a linear 
system of equations, thus making the calculations slow and expensive). The spline solution. 
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however, is technically very good, and if the cost/benefit is worth, maybe in certain applications 
one should really use the splines. 

The solution proposed in this work is the weighted splines method (Costa, 1980; Pereira and 
Vijaykumar, 1984) that has the advantage of low computer cost since it docs not have to solve a 
linear system of equations, produces a 3rd degree polynomial based on four data points calculated 
through a single formula and is very convenient to be computer implemented. As a further advantage, 
it is possible to impose the values of the first derivatives, in or^er to account for sharp 
boundaries, without increasing the computational effort or to keep the first and second derivatives 
of the solution function continuous as determined by the data poirts (Pereira and Vijaykumar, 1984). 
In the examples no sharp boundaries were considered. 

The weighted spline method is described and a few examples ai'e presented to illustrate it. 

The advantages and shortcomings are discussed in the conclusions. 

2. WEIGHTED SPLINES METHOD 

The interpolation .n each interval i is obtained using only four points. Given n pair of 
points (x^, y^l such that i=1, 2, 3, ..., n and x. < the following hypotheses are assumed: 

(i) f(x.) » y.; 


(ii) f'(x) and f"(x) are continuous; 


(iii) in each interval [x^, x^^^] f(x) mu:t be defined by a polynomial wliose 
degree does not depend on n; 

(iv) in each interval [x^, x^^^], f(x) must be defined only as a function of 

X. and y., where j=i-1, i, i+1, i+2. 

J J I 

Initially, two second degree functions hj and hj are defined such that hi interpolates by 
*i-i, X. and x. and h» interpolates by x., x. and x. . The following conditions are given to 
these functions: 


hi(x._j) - y^_j-, hi(x^) - y.; hi(x.^^) 

hj(x.) = y.; h,(x.^j) = y^^^; hi(x.^J 

Now consider a function f(x) defined as 

f(x) = Pi(x) hi(x) ♦ Pi(x) hj(x), 
where Pi(x) and Pj(x) are weight functions. 



( 2 . 1 ) 

( 2 . 2 ) 


(2.3) 


For x-x^ and x=x^^j the function f(x) in Equation 2.3 yields, respectively: 
f(x^) * (Pj(x^) + Pi(x^)) y^. 


• ^(*i+i^ ■= ('’‘(’‘i+x^ * P^(’'i+i)) ^i+x- 

From Equations 2.1 and 2.2 it is known that f(x.) c y,. With this condition the following can be 
obtained: 

Fi(x^) + Px(x,.) * 1 ; 


t 
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Deriving Equation 2.3, one may obtain: 

f'(x) - Pi(x) hx(x) + Pi(x) h|(x) + Pi(x) h 2 (x) + P 2 (x) h^x). (2.4) 

In order to guarantee the continuity of the irst derivative, the following has to be 


satisfied: 

f'(x.) c hl(x.). 

The following conditions satisfy the above: 

Pl(x.) + Pi(x.) = 0, (2.5) 

Pi(x^) - 1, (2.6) 

P2(x.) = 0, (2.7) 

Pl(x^^,) + Pi(x^^^) - 0, (2.8) 

Pi(x^^j) « 0, (2.9) 

P2(x^^,) = 1. (2.10) 


For f"(x), the second derivative of f(x), the continuity must also be guaranteed and the' 


following conditions are necessary: 

Pi*(x^) ♦ Pr(x^) - 0, (2.11) 
Pl(x^) =0, • (2.12) 
Pi(x.) = 0, (2.13) 
P;’(x^^j) + Pi'(x^^j) = 0. (2.14) 
PJ(x.^^) . 0. (2.15) 
Pi(x.^J = 0. (2.16) 


Changing our interval to [0,1], and changing the variable from x tc z, t'.e following 
relations hold for x t [x^, x^^^] and for z c [0, 1]: 


*i ‘ *i*i 

The above system of equations in z becomes: 
Pi(0) = 1; Pj(0) . 0. 

Pi(1) - 0; P^d) . 1. 

P;(0) = 0; Pi(0) . 0. 


PI(1) . 0; Pl(1) - 0. 
Pl'(O) + Pr(0) - 0. 
pi'd) + Pi'd) - 0. 
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Pi(z) and p 2 (z) can be expressed as fourth degree polynomials! 


(2.i7) 


Pi(z) = a^z** + hjZ* + CjZ^ + djZ + e^i 

Pj(z) = a 2 z" + b 2 z’ + C 2 Z^ + d 2 Z + ea. (2.18) 

Substituting Equation 2.18 in Equations 2.17 and solving the system, the following is 
obtained. 

Pi(z) • (1 + |)z‘* - tz’ + (| - 2)z* + 1, 

P 2 (z) . -d +|)z- + tz’ - - 2)z^ 

where t is a parameter to be determined. It can be noticed that 
Pi(z) - 1 - P 2 (z) or PjCz) = 1 - Pi(z). 

The Expression 2.3 can be written: 

f(z) = Pi(z)hi(z) + (1 - Pi(z))h 2 (z), 

which is the interpolating polynomial. It was denwnstrated by Costa (1980) that if t c [-8, 4) the 
oscillations are avoided. 

Returning from z to x, the final formula is 

f(x) * Pj(x)hi(x) + P 2 (x)h 2 (x), 


where 





P2(X) = 1 - Pl(x) 
t c [-8. 4). 

One can also notice that when t « -2, the above functions are reduced to third degree 
polynomials. 

The functions hjx) and h 2 (x) are: 


h,(x) c 


^’‘i.i - "i^ (’‘i-x * *i.J 


(X - X, .i) (X - X^^i) y ^ 

(x. . X._J (X. - 


^i-l ^ 
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^ (x - Xj.,) ( X - x^) 


hj(x) 


(x - xui) (x - xi-;,) ■ , (x - Xj) ( X - xuJ ^ 

(x^ - X^^J (x^ - X^^,) ' (x^^^ . x^) (x.^^ . x^^^) 


^ (x - Xi) (x - X^^i) 


yi«- 


3. CONCLUSIONS 


In order to test the method, a computer program was developed having as input a pixel matrix 
8x8, with each pixel in the range 1(one) to 255 (two hundred and fifty-five). The level 0 (zero) 
was considered noise or interference. 

Depending on the target, the weighted splines method can produce better or worse results. If 
the missing line is the upper or lower border, the averages method does not give results and the 
weighted splines method docs. If there is more than one adjacent lines lost, again the averages 
method does not work, and the weighted spline, method docs. 

To cti,)Osr the better reconstruction, the elements of the vectors for the reconstructed lines 
were simply a^-ded, the lower the sum, the better the reconstruction. 


Two examples were considered: 

1) Using the artificial image below (pixel values from 1 to 255), 
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the lines 3(thrce) and 7(scven) were removed and substituted for zero (assumed to be noise), and 
then the average method was applied for reconstruction. The results v/ere: 

for line 3: 75 75 75 75 75 75 75 75 

for line 7: 77 67 162 140 171 167 81 41, 

and the sum of the lino 3 produced 200, while the sum of the line 7 produced 441. 

Applying the weighted splines method to the example, it was possible to obtain: 

for line 3: 44 44 44 44 44 44 44 44 

for line 7: 61 21 144 152 162 241 63 37, 

and the sums for lines 3 and 7 resulted in 48 and 408, respectively. 


In this case the weighted splines method showed to be better. It has to be said that line 3 
was chosen such that average method did not work. However, for line 7 a random generation v/as used 
to produce the pixel values and, as one can notice, the weighted splines method v/a> better. 


2) Using a portion of a Thematic Mapper Image (this image belongs to Southern Bahia (Brazil) region 
with oscillating pixel values). 


- 6 - 


64 

71 

64 

74 

59 

51 

59 

69 

59 

59 

59 

56 

59 

48 

51 

69 

59 

56 

61 

59 

61 

61 

59 

71 

61 

61 

71 

64 

61 

59 

69 

71 

69 

59 

64 

64 

66 

61 

59 

59 

94 

94 

94 

66 

74 

74 

71 • 

61 

145 

94 

74 

59 

87 

105 

107 

102 

59 

48 

38 

54 

87 

102 

107 

105 


the lines 3 and 7 were removed and substituted for zero. The following results were obtained using 
the averages method and weighted splines method: 


Averages method: 


line 3: 60 
line 7: 76 
sum of line 3: 
sum of line 7: 


60 65 

71 66 

21 
159 


Weighted Splines Method: 


line 3: 57 
line 7: 90 
sum of line 


58 65 

90 85 

: 24 

sum of line 7: 143 
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The averages method presented a better result for line 3, and the weighted spline method 
presented a better result for line 7. 


Further tests have to be made for better and precise conclusions. However, it is expected 
that in near-random fluctuations of pixel variation, the two methods v/ould yield similar results. 
In nonlinear variations (as in the case of most number of images), the weighted splines method 
would be better. In linear variations the averages rcthod presents better results. 

As far as computer effort is concerned, both the averages method and the weighted splines 
method present the same complexity (order n), although the former is slightly faster. 
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